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The density matrix of a nonrelativistic wave-packet in an 
arbitrary, one-dimensional and time-dependent potential can 
be reconstructed by measuring hydrodynamical moments of 
the Wigner distribution. An n-th order Taylor polynomial in 
the off-diagonal variable is obtained by measuring the proba- 
bility distribution at n + 1 discrete time values. 
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This Letter presents a new and general method for re- 
constructing the density matrix of a massive particle in 
an arbitrary, one-dimensional and time-dependent poten- 
tial. Such a general method may seem called for, e.g., in 
the reconstruction of the quantum state of particles in 
anharmonic, time-dependent Paul traps. The method is 
based upon measuring the position probability distribu- 
tion for a discrete number of time values in a short time 
interval. Surprisingly, the method follows almost imme- 
diately from known results. 

The diagonal of the density matrix can be retrieved by 
observing a single probability distribution. The state re- 
construction problem essentially consists in obtaining the 
offdiagonal elements. Decoherence and the approach to 
the classical regime is characterized by a vanishing of the 
offdiagonal elements [Q. The method presented here is 
constructed so that the density matrix is retrieved for in- 
creasing values of the off-diagonal variable by increasing 
the number of discrete time values for which the proba- 
bility density is observed. 

It was shown by Madelung Q that quantum mechan- 
ics can be reformulated in a form resembling a hydro- 
dynamical description. He reformulated the Schrodinger 
equation as two coupled and nonlinear equations for the 
"hydrodynamical" moments of probability density and 
probability current density. Whereas these two moments 
can describe a pure state ||, the situation is more com- 
plicated for mixed states. 

A somewhat analogous situation is found in classical 
statistical mechanics. Hilbert [0] demonstrated that the 
phase space distribution for a system in local thermal 
equilibrium can be expressed as a functional of the den- 
sity, the current density and the kinetic energy density. 
These are the three lowest order velocity moments of the 
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phase space distribution. Thus, in thermal equilibrium 
the phase space distribution is equivalent to it's three 
lowest order moments. This situation has sometimes 
been called the Hilbert paradox ||. In general, though, 
an infinite set of velocity moments is equivalent to the 
full phase space distribution. These moments are cou- 
pled through an infinite set of differential equations. In 
the case of thermal equilibrium this set is truncated, and 
a finite and closed set of equations is obtained. 

The similarity between quantum mechanics and sta- 
tistical mechanics goes beyond the observation made b 
Madelung, which is restricted to pure states. Wigner 
showed that quantum mechanics can be reformulated in 
terms of a quasi phase space distribution. This distribu- 
tion, which is the Fourier transform of the density matrix, 
shares most of the properties of a classical phase space 
distribution. One important exception is that it may take 
on negative values. It can be shown j?J that the proba- 
bility density and the probability current density are the 
first two velocity moments of the Wigner distribution. 
This is in complete analogy with classical statistical me- 
chanics. An infinite hierarchy of velocity moments can 
be derived from the Wigner distribution, and they are in- 
terconnected through an infinite set of coupled equations, 
much like in classical statistical mechanics. 

However, quantum mechanics can offer more. It has 
been shown p|-|ll| that when the density matrix in posi- 
tion representation is expanded as a Taylor series in the 
offdiagonal variable, the coefficients of this expansion are 
velocity moments of the Wigner distribution. This in fact 
solves the moment problem both for classical statistical 
mechanics and for quantum mechanics. With "the mo- 
ment problem" we here mean the problem of expressing 
a distribution in terms of it's moments. This is a classical 
problem in statistics, and it was first raised in the con- 
text of quantum mechanics by Moyal It was further 
explored in Ref. Jl2| . Recently, the density matrix was 
expressed in terms of normally ordered moments []l3f . It 
has also been shown Jl4| that the normally ordered mo- 
ments can be calculated from the measured quadrature 
distribution. Since convergent expansions may be found 
jyl, this gives a method for reconstructing the state of 
a radiation field or a particle in a harmonic oscillator 
potential. 

The possibility of measuring quantum states has at- 
tracted a lot of attention in recent years. In particular, 
the method of homodyne tomography |L6 17 has con- 
tributed to this interest. It has been used to reconstruct 
the quantum state of radiation fields |l8| as well as mate- 
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rial particles propagating in free space fll9[ . In homodyne 
tomography the state is retrieved from a parameterized 
probability distribution. Ideally, a continuous range of 
parameter values should be used. In the case of opti- 
cal homodyne tomography, this parameter is the value 
of a reference phase [Il7j, whereas for material particles 
it might be a time value . In another recently de- 

veloped reconstruction scheme, a method for the direct 
probing of the Wigner distribution has been found [p2| . 
For a specific parameter value (in this case, the amplitude 
and phase of a probe field) a certain region of phase space 
is retrieved. This method has been used to reconstruct 
the first negative Wigner distribution ^3|. Numerous 
other reconstruction methods have also been found 0| . 

Recently, it has been shown that the density matrix of 
a wave-packet in an arbitrary one-dimensional potential 
can be reconstructed by observing the time-evolution of 
the position probability density |pl| , p5[ . In these meth- 
ods, the eigenstates of the Schrodinger equation are first 
found for the potential in question. The position proba- 
bility density should ideally be observed over an infinite 
time interval, although methods have been considered for 



obtaining a finite observation time 25 26y. 

The density matrix is the Fourier-transform of the 
Wigner distribution || 

/oo 
dpe 2 ^/ n W(x,p,t). (1) 
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where /„ are "hydrodynamical" moments of the Wigner 
distribution [HI 
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Since the Wigner distribution is real, these moments are 
also real, fo trivially is the probability distribution in 
position representation, whereas fi/m is the probability 
current density fij. Generally it is not possible to give 
these moments a classical hydrodynamical interpretation. 
Thus, e.g., the moment fi may take on negative values 
for certain negative Wigner distributions |^7j . 

It follows that the unique Taylor expansion of the den- 
sity matrix in the off-diagonal variable y is |)] [lTj 
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We see that the real part contains only moments /„ of 
even order n, whereas the imaginary part contains only 
moments of odd order. 

Clearly, if we are able to measure the moments /„, 
we have a state reconstruction scheme. To this end, we 
recall the equation of motion of the Wigner distribution 
for a particle with mass p in a one-dimensional, time- 
dependent potential V(x,t), PjO] 
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We multiply both sides with p n and integrate over all 
momentum space. In this way, we obtain the infinite set 
of coupled equations pX| , pH 
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For n — we retrieve the well known conservation 
equation for probability. By integrating this conserva- 
tion equation, the probability current density can be ex- 
pressed in terms of the time-derivative of the cumulative 
position probability f28f| (for one-dimensional systems, 
that is). The idea of the present reconstruction method 
is simply to generalize this procedure to arbitrary mo- 
ments. This will yield an iterative scheme. We therefore 
integrate Eq. (Q) over the position variable and obtain 
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We may divide this expansion into a real and an imagi- 
nary part by 



The moment f n +i is expressed in terms of lower order 
moments only. Therefore an arbitrary moment can be 
recursively calculated from the zeroth order moment, the 
probability density. This recursion relation gives an al- 
gorithm for reconstructing the density matrix. The al- 
gorithm can be used directly on the experimental data. 
In order to find /„, we must know the time derivative 
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of fn-i- Therefore, f n -i must be observed for at least 
two different time values. This again requires that / n _2 
is known for three time values. Recursively, it follows 
that /„ can be found by measuring fo for at least n + 1 
different time values. 

To illustrate the convergence of the Taylor series (|4|), 
consider the unnormalized superposition state 
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The corresponding density matrix is 



(x + y\p\x -y) = 2exp 



x 2 + y 2 



2a 2 

x [cos(2k x) + cos(2fc y) 
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Note that this density matrix is real. This means, ac- 
cording to Eq. (||), that f n vanishes for all odd n. In 
Fig. [I] the Taylor polynomial 
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has been plotted for different orders N for the parame- 
ter choice er = 1/V2 and k = 2\[2. The highest order 
polynomial is p^g (Fig. [l] c), which involves moments up 
to /36- It would be obtained by measuring the position 
probability distribution fo for 37 different time values us- 
ing perfect detectors. It differs negligibly from the exact 
density matrix within the chosen plotting region. 

As we can see from Fig. [lj one in effect probes the den- 
sity matrix further away from the diagonal by increasing 
the number of discrete time values for which the prob- 
ability distribution is observed. If, after retrieving the 
Taylor polynomial (|ll]) to a certain order, one finds that 
the density matrix goes to zero, one may use an addi- 
tional number of measurements at other time values to 
check the consistency of the data. 

The classical limit is often associated with taking Ti — > 
0. In this case the equation of motion (||) reduces to a 
classical Liouville equation. But for the purpose of re- 
constructing the Taylor polynomial (0), it is vital that 
h should be considered finite. Otherwise, assuming finite 
hydrodynamical moments /„, every term in the polyno- 
mial of order higher than zero diverges. The specific nu- 
merical value of h is of no importance in this respect, 
since a change of fi only implies a rescaling of the off- 
diagonal variable. 

The principles outlined in this Letter can also be em- 
ployed to other expansions of the density matrix. Start- 
ing with a density matrix in the momentum representa- 
tion, we might have expanded it in terms of the moments 
dxx n W(x,p, t). However, the corresponding set of 
recursion relations is generally more complicated in this 
case. Other expansions of the density matrix, which con- 
verge more rapidly for nearly classical states EE, might 



also be used for state reconstruction using similar tech- 
niques. It is also easily adapted in other areas such as 
the reconstruction of quantum optical states. 
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FIG. 1. The real part of the Taylor polynomial pjv of the 
density matrix for a) N = 10 b) N = 20 and c) N = 36. 
It is seen that the density matrix is retrieved for increasing 
values of the offdiagonal variable as the number of discrete 
time values are increased. p36 is almost indistinguishable from 
the exact density matrix within the chosen plotting region. 
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The expansion (||) of the density matrix in terms of 
quasi-hydrodynamical moments can be generalized to 
systems with a higher number of dimensions. However, 
it may not be straightforward to generalize the recursion 
algorithm (||). This is basically due to the fact that the 
current density is not uniquely determined from the time 
derivative of the probability density for systems with two 
or more dimensions. 

In conclusion, a method was found for reconstructing 
the density matrix of a particle in an arbitrary, time- 
dependent potential. The method was based upon a Tay- 
lor expansion of the density matrix in the off-diagonal 
variable. The coefficients /„ in this expansion are hydro- 
dynamical moments of the Wigner distribution. A re- 
cursive algorithm was found for calculating an arbitrary 
moment /„ from the zeroth order moment, the probabil- 
ity distribution. In general, an n-th order Taylor poly- 
nomial of the density matrix can be found by observing 
the probability distribution at n+ 1 discrete time values. 
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